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Following the observational evidence for cosmic acceleration which may exclude a possibility for the uni- 
verse to recollapse to a second singularity, we review alternative scenarios of its future evolution. Although 
the de Sitter asymptotic state is still an option, some other asymptotic states which allow new types of sin- 
gularities such as Big-Rip (due to a phantom matter) and sudden future singularities are also admissible 
and are reviewed in detail. The reality of these singularities which comes from the relation to observational 
characteristics of the universe expansion are also revealed and widely discussed. 
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1 Introduction 

It is generally agreed that we have now enough evidence for the past Hot-Big-Bang universe. Its main 
observational support relies on the following facts: 

• The universe expands, i.e., all the galaxies move away from each other according to the Hubble law 
1 1 1 and they experience cosmological redshift z according to the relation: 

i+z = 454, a) 

a(t e ) 

where a(to) is the scale factor at the time of observation, while a(t e ) is the scale factor at the time of 
emission of light by a galaxy. 

• Element abundance in the universe is: hydrogen 75%, helium 24%, and other elements 1%. In partic- 
ular, the amount of helium is larger than it is possible to be produced in stars, and the only solution to 
this problem is to assume that its abundance is primordial \ 2\. 
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• Cosmic Microwave Background (CMB) - photons once were in thermal equilibrium with charges 
which further decoupled and formed thermal background with blackbody radiation spectrum with 
temperature T — 2.7 K Q. The information about the density fluctuations Sp at the decoupling 
epoch was imprinted in the temperature fluctuations according to the formula 

cx ^oclCT 5 . (2) 
T p 

However, as well as about the past, it is interesting to ask about the future of the universe. The questions 
which naturally arise are as follows: 

• What type of future evolution will we generally face? 

• Will the universe expand forever? Will it expand fast, faster, slower . . . ? 

• Will we face any dramatic change of our future evolution? 

• Is it likely that we face an unexpected end of our future evolution? 

• Is there a barotropic equation of state p(t) = wg(t), w = const, (p- the pressure, g- the energy 
density) valid throughout the whole evolution, or, perhaps w — w(t) so that the pressure can be 
expanded in series as 

dp 1 d^p 

p = Po + -j- lo {q- go) + 2\~d^ 1° ( e ~ Q ^ 2 + £°) 3 ] 

, Po, \ 1 Po0o ~ Po0o , \2 , n\i ,cn 

= Po + — [Q- Qo ) + 7^ ^ [8- 0o) +0 Ug- g Q , (3) 

go 2! g 6 

where index "0" refers to a quantity taken at the current moment of the evolution t = to- 

• To what extend we are able to determine wq = Po/ go, W — (dp/dg) |o, ■ ■ ■ ? 

It is generally believed that these questions may, in general, be addressed within a more fundamental 
framework than Einstein's general relativity theory, i.e., within the framework of fundamental theories of 
all physical interactions such as superstring, brane and M-theory ffl|5|. 

2 Empty future, Big- Crunch and phantom-driven Big-Rip 

We start our discussion from the Einstein's field equations for the homogeneous and isotropic Friedmann 
universe in the form (we have assumed that 8ttG = c = 1) 

0(t) - 3 ( ^ + *) , (4) 
( a a 2 K\ 

P(t) = - 2- + - + - , (5) 
\ a a z a z ) 

where a(t) is the scale factor, K = 0, ±1 is the curvature index. These two equations contain three 
unknown functions a,p, g. In order to solve the system one usually assumes the equation of state of a 
barotropic type, i.e., 

p(t) = wg(t) (6) 

with w = const, which leads to the three solutions - each of them starts with Big-Bang singularity in which 
a — > 0, g,p — * oo - but only one of them (of K = +1) terminates at the second singularity (Big-Crunch) 
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where a — > 0, g,p — > oo while the other two (K = 0,-1) continue to an asymptotic emptiness g,p — > 
for a — > oo. Besides, at least one singularity (e.g. Big-Bang) appears provided the strong energy conditions 
of Hawking and Penrose |6) 

R^V^V" > 0, V - a timelike vector , (7) 
- Ricci tensor) is fulfilled. In terms of the energy density and pressure it is equivalent to 

Q + 'ip > 0, Q + P > . (8) 
From and (|3J one has 

a 4nG , 

- = ^(Q + 3p), (9) 

which together with (|8j means that 

a > 0, (10) 

so that the universe decelerates its expansion. 

However, the observations of type la supernovae |7| in 1998 gave the evidence for 

S<0, (11) 

which means that the universe currently accelerates its expansion. Taking into account l|9} it means that 
there exists negative pressure matter (dark energy, quintessence) 

P<~\q (12) 

in the universe which drives this acceleration. More precise fit to the data shows that at least 70% of matter 
in the universe now has negative pressure and it is neither the visible matter, nor the dark matter, which 
can sum up to the remaining 30%. Some possible sources of negative pressure had already been suggested 
before 1998 and they were: the cosmological constant w = —1 |8|, cosmic strings w = —1/3, and domain 

waiis w = -2/3 sua. 

In fact, this result gives more evidence for an asymptotic emptiness in the future of the universe de- 
scribed by the de Sitter model in which q — > when a = exp [(-y/A/3)t] — > oo. From \\2\> it is obvious 
that the strong energy condition is violated. 

Last but not least, the most recent observations of type la supernovae show II 1111 21 that the pressure may 
not only be negative, but it may also be very strongly negative. In other words, there is not an observational 
barrier onto the amount of negative pressure and it is also very likely that 

p < —q, or w < —1 . (13) 



This is called phantom \ 13 14 1 and it is the dark energy of a very large negative pressure which violates 
all the remained energy conditions, i.e., the null 

T^Wk" > 0, W - a null vector , i.e., Q + p>0, (14) 
(T^ v - an energy-momentum tensor), the weak 

T^V^V" > 0, V 1 - a timelike vector , i.e., Q + p>0, p>0, (15) 
and the dominant energy 

T lu/ V lt V v > 0, T^V^ - not spacelike , i.e., \p\<Q, Q>0. (16) 

The violation of these energy conditions seems to be the most difficult problem for the models admitting 
w < — 1 because of the emergence of both classical and quantum instabilities 1 15 1 and there are some 
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suggestions it should be avoided 1 16 1. Despite violation of all the energy conditions phantom allows for a 
new type of future singularity - Big-Rip which characterizes by a blow-up of both the energy density and 
the pressure (g,p — > oo) together with the infinite size of the universe a — > oo. 

In order to make this divergence clear, let us briefly study the basic properties of phantom cosmology. 
From the conservation law we have 



Q CX 



-3(«.+i) for w > _1 (17) 



and 



g cx a?l w+1 \ for w < -1 (phantom). (18) 

From Eq. Jl 8i it is clear that the singularities appear at infinite values of the scale factor for phantom 
cosmologies. This is also obvious after studying the simplest solution of J4} for one general fluid p — wg 
only, which gives 1141 

a(t) = a | t |3T^ftt ; (19) 
for ordinary (w > — 1) matter, and 

o(t) - a 1 1 |"W , (20) 
for phantom w + 1 = — | w + 1 |<0 and 

q cx r 2 . (21) 

In other words, taking w = —4/3 (phantom) one has a(t) — > oo and g — ► oo if < — * 0, while a(t) — > and 
g — * if £ — * oo. On the other hand, in a standard case w = —2/3, for example, one has a(t) — > and 
g — > oo if i — * 0, while a(i) — * oo and p — > if t — > oo. 

It is worth noticing that both the non-phantom matter (— g < p) and the phantom matter (p < —g) may 
be mimicked by a scalar field <j> with some potential V(<p) with the effective energy density and pressure 

q = ±^ 2 + V(0), (22) 

p = ±^ 2 + V(0), (23) 

where the plus sign refers to the non-phantom matter and the minus sign refers to the phantom. From the 
formulas (I22M 23>. it follows that phantom can be interpreted as a scalar field with negative kinetic energy 
(a ghost). 

Another interesting remark can be extracted from the Eqs. @-(|5j and dl7t -( ll8t if we admit shear 
anisotropy a^/a 6 (<7 = const.) and consider nonisotropic Bianchi type IX models. Namely, for w < —1, 
the shear anisotropy cannot dominate over the phantom matter on the approach to a singularity when 
a — > oo, i.e., we have 

2 

ga 3\w+i\ > 1SL f or a ^oo (24) 

and this prevents the appearance of chaotic behaviour of the phantom cosmologies of the Bianchi type IX 

Bearing in mind the fact that the Big-Bang/Big-Crunch singularity appears for a — > while the Big-Rip 
singularity for a — > oo one may suspect a kind of duality between the standard matter (p > 0)/quintessence 
(— g < p < 0) models and phantom (p < —g) models which is present in the low-energy-effective 
superstring theory 111 81 . Indeed, there is such a duality, called phantom duality, which explicitly reads as 
fit) 

w + l <-> -(tu + 1), a<->~. (25) 

a 
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This duality can easily be seen if we rewrite the system of equations @-(|5j in the form of the nonlinear 
oscillator 

D 2 

X - —AX + D(D - \)kX x - 2 l D = 0, (26) 
after introducing the variables 

X = a D{w \ D(w) = ^(l + w). (27) 

It is obvious to notice that Eq. (I26> preserves its form under the change 

D <-> —D . (28) 
Alternatively, this invariance takes form (H = a/a) llT9l 

H ^ -H , g + P^-(g + p)- (29) 

In fact, there is a richer symmetry of the field equations which includes brane models called phantom 
triality Eol . 

The simplest way to consider these dualities is to look at the solutions J17I - (I18I >. In both cases there is 
a curvature singularity at t = 0, but in the former case it is of a Big-Bang type, while in the latter case it 
is of a Big-Rip type. From the observational point of view it is reasonable to choose the solution dl9t for 
positive times t > 0, and the solution J20i for negative times t < 0. Another example of a phantom model 
with an explicit phantom duality is (for flat models with walls w = —2/3, phantom w = —4/3, and A < 

El) 

a w = ^sin^=HA|t) , (30) 

where 




a ph = S mL^|A|t] , (31) 



D w = 1/2 = -D ph , (32) 
so that we have 

a w = %h ■ (33) 

It is obvious that the evolution of a w begins with Big-Bang and terminates at Big-Crunch while the evolu- 
tion of a p h begins with Big-Rip and terminates at Big-Rip (cf. Fig|2jl. 

However, from the point of view of the observations (which support Hot Big-Bang models), the most 
interesting models are "hybrid" models which begin with Big-Bang singularity and terminate at Big-Rip. 
Both types of matter (standard and phantom) are present during the whole evolution of the universe but 
an early evolution is dominated by the standard matter, while phantom dominates the late evolution of 
it. Of course, this means that there must have been a change during the evolution from deceleration to 
acceleration and this also might have happened quite recently (we will come to this point later). An explicit 
example of such an evolution for the dust (w = 0) and phantom (w = —4/3) model is given in terms of 
Weierstrass elliptic functions as 



a\ iQ m0 V( V ) Vo V; (34) 



a J AQ, p hoV 2 {-q) - fi m0 

where (77) is the Weierstrass elliptic function, H, m o,n p ho the density parameters of dust and phantom 
respectively. 
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Fig. 1 The typical phantom solution which begins at an initial Big-Rip singularity (a — » oo, g, p —> oo) and terminates 
at final Big-Rip singularity (a — > oo, g,p — > oo). Similarly as in the case of Big-Bang-to-Big-Crunch evolution many 
cycles are presumably admitted. 



a 




t 



3 Sudden future singularities of pressure and generalized sudden future 
singularities 

Big-Rip may appear in some future time t = tsR analogously to Big-Crunch (which may appear in some 
future time t = tBc), but because of growing acceleration it is sometimes called "sudden". However, we 
may have something more exotic in the future evolution of the universe - a singularity which presumably 
appears quite unexpectedly and does not violate all the energy conditions. The hint which allows for such 
a singularity is that we release the assumption about the imposition of the equation of state, i.e., we do not 
constrain pressure and the energy density in @ and @ by any equation like the one in ©. This enables 
quite independent time evolution of these physical quantities. 
Suppose that we first choose the form of the scale factor as 1211 

o(t) = A+(a s -A)(l) -a(i-I) , (35) 



(A = const., a s = a(t s )) with its time derivatives 

m , . , / t \ m . n ( t \ n 1 

- " u { -- A Au) <7 ( ~ ~J • ^ (36) 

m (to — 1) . . / t \ n(n — 1) / t x 



Choosing 

1< n < 2, < to < 1 . (38) 

we notice that the scale factor J35i vanishes and its derivatives d36l >- d37l > diverge at t = leading to a 
divergence of g and p in (|4j-(|5j (Big-Bang singularity). On the other hand, the scale factor J35i and its first 
derivative d36l > remain constant while its second derivative d37l > diverge leading to a divergence of pressure 
in @ only, with finite energy density ©, i.e., 

a = const., d = const., d — > — oo, (39) 
q = const., p — > oo, for t — > t s . (40) 

We then conclude that the future singularity appears and since it does not show up in the evolution of 
the scale factor, it is called "sudden" future singularity (SFS) |21 1. This singularity violates the dominant 
energy condition d!6t \22\, but does not violate any other energy condition. 
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Sudden future singularities may also be present in the inhomogeneous models of the universe where 
they may violate all the energy conditions. An explicit example of such a situation was given |23| for the 
inhomogeneous Stephani universe [24 1 with no spacetime symmetries with metric 



, , a 2 a 2 

a z V z 



2 



dt 2 + — [dx 2 + dy 2 + dz 2 ] , (41) 



where 



V(t,x,y,z) = l + h(t){[x-x Q (t)} 2 + [y-ya(t)} 2 + [z-z Q (t)} 2 } , (42) 



and a, xq, yo, zq are arbitrary functions of time. It emerges that sudden future singularities are similar to 
finite-density singularities of pressure (FD) which appear in these Stephani models since they also do not 
admit a priori any equation of state. Besides, sudden future singularities are temporal pressure singularities 
which means that they may appear at an instant t = t s of the evolution, while finite-density singularities 
are spatial pressure singularities and they may exist somewhere in the universe nowadays. Sudden future 
singularities were also proven to exist in anisotropic models 11261 . 
If we consider a general time derivative of an order r for (I35> . i.e., 



(r) _ m(m — l)(m — 2)...(m — r + 1) . . / t \ rn 

v. {as ~ 'VsJ 



(43) 

v s \ [ s J 

and replace the condition d38i by 

N < n < N + 1, < m < 1 , (44) 

we realize that for any integer N + 1 = r we have a singularity in the scale factor derivative (cf. j43». 
and consequently in the appropriate pressure derivative p( r ~ 2 '. This, for any r > 3, gives a sudden future 
singularity which fulfills all the energy conditions including the dominant one 1251 1221 . These singularities 
are, for example, possible in the theories with higher-order curvature quantum corrections 1 27 1. 

In fact, sudden future singularities are determined by a blow-up of the Riemann tensor and its derivatives 
rather than the scale factor (or its inverse) itself. This has an important consequence on the nature of these 
singularities in comparison with Big-Bang or Big-Rip singularities. Namely, the geodesies do not feel the 
sudden future singularities at all, since geodesic equations 

dt Y a P 2 + KL 2 

— J = A + — . < 45 ) 

(46) 

(47) 



dr J a(t) 
dr Picoscf) + P 2 sin i 



dr o(t)/(r) 
d(f) L 
dr a{t)r 2 ' 



are not singular for a s = a(t s ) = const. 1288 . Here (t, r,0 = 7r/2, <f>) are coordinates of the Friedmann 
universe, r is the proper time, K = 0, ±1, A, L, P 2 = P 2 + P| are constants, and f 2 (r) = 1/(1 — Kr 2 ). 
On the other hand, the geodesic deviation equation 

+ R% 6 u^n?u s = , (48) 

where n a is the deviation vector, and u 13 is the four-velocity vector (A - an affine parameter), feels this 
singularity since at t = t s we have the Riemann tensor R a p 7 g — > oo. This means point particles do not see 
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sudden future singularities, but extended objects do in the sense that they may suffer infinite tidal forces 
although they may not be torn on crossing these singularities. 

Besides, the spacetimes with sudden future singularities do not lead to geodesic incompletness - they 
are weak singularities according to the definitions of Tipler |29| and Krolak |30|. In other words, sudden 
future singularities are not final state of the universe 1281 . This situation is somewhat similar to what 
happens with finite density singularities in inhomogeneous Stephani models 1 3 1 1 although there instead of 
the extension of geodesies one has the extension of the hypersurfaces of constant time throughout these 
singularities. 

In order to conclude, one has to stress that sudden future singularities and finite density singularities are 
totally different from Big-Rip singularities due to phantom since Big-Rip singularities are felt by geodesic 
equations and the evolution cannot be extended behind them. In other words, Big-Rip (like Big-Crunch) is 
really the final state of the evolution while sudden future singularity is not. 

4 Statefinders and the diagnosis of the future state of the universe 

From the discussion of the previous sections it is clear that all the types of singularities (Big-Bang, Big- 
Crunch, Big-Rip, sudden future singularities and generalized sudden future singularities) are related to a 
blow-up of the scale factor and its time derivatives 

a(t),a(t),a(t),a{t)... , (49) 

and so to a blow-up of the energy density and pressure derivatives 

Q{t),Q{t),...,p(t),p(t),... . (50) 

Formally, we may translate all that into the observational characteristics of the expansion (statefinders) 
such as: the well-known Hubble parameter 

H = °, (51) 
a 

and the deceleration parameter 

9 = -7777 = ~77 . ( 52 ) 



(53) 



k = - — - = -—> ( 54 ) 



together with the new characteristics: the jerk parameter J32 

1 a a a 2 
j = Fa = ~~dF ' 
and the "kerk" (snap) parameter I33II34I 

1 o 

the "lerk" parameter 

1 = -m— = —^i (55) 

and "merk", "nerk", "oerk", "perk" etc. parameters, of which a general term may be expressed as 

* W - (-D i+1 4— = (-l) m ^^> (56) 



and its time derivative reads as 

(x®y = H U(q + l)a;W _ LW) + X (A\ 



(57) 
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Future predicted ("sudden") blow-up of statefinders may easily be linked to an emergence of future 
singularities. In particular, this may be the proper signals for sudden future and generalized sudden future 
singularities. 

The blow-up of statefinders can be read-off redshift-magnitude relation up to an appropriate order in 
redshift z after using the redshift-magnitude relation applied to supernovae data 1 3 31 1341 



m 



M = 5 log 10 {as) 



5 log 10 H [ - log 10 e 



{l-q )z+- 



TT 77<?0 



— jo — J — 



(58) 



1 



[2 jo (8(70 + 5) - 2k - q {7ql + llgo + 23) + 25 



24 

+4fW2 go + l)]z 3 + 0(z 4 )} ■ 
where m is the observed magnitude, M is the absolute magnitude, c is the velocity of light, and Qko = 
Kj (Hy ag). From J58l > it is clear that the jerk appears in the second order of the expansion and the "kerk" 
appears in the third order of this expansion. In particular, the simple signal for Big-Bang, Big-Crunch 
or Big-Rip would be | H Q |— > oo; for SFS would be | q |— > oo, while for generalized SFS would be 
| j |-> oo, | k |-* oo etc. 

As it was mentioned in the Introduction, statefinders may be useful in determination of the current status 
of an equation of state Q for the matter in the universe. 

The most recent analysis of supernovae type la data 1 12 1 shows that 

jo > 0, (59) 

giving no hint about ko. Apart from that it shows that acceleration of the universe started quite recently at 
redshift z = 0.46 ±0.13. 



5 SFS avoidance - generalized energy conditions 

From the analysis of the previous Section|4]we learned that Big-Crunch singularity in future may emerge 
despite all the energy conditions are fulfilled. On the other hand, Big-Rip may emerge despite the energy 
conditions are not fulfilled. Further, SFS emerges when only the dominant energy condition is violated, 
while generalized SFS do not lead to any violation of the energy conditions. This means that an applica- 
bility of the standard energy conditions of Hawking and Penrose 1 6 1 to cosmological models with more 
exotic properties is limited and so not very useful. 

Because of that, one may consider some more general or different energy conditions which may be 
helpful in classification of singularities in contemporary cosmology based on the fundamental physical 



theories such as superstring theory, brane theory, or M-theory 1 36 37 38 39 40l l41l . 
First, let us consider the following higher-order dominant energy conditions 1 35 1: 

Q > 0, ~P<g<p, (60) 
Q > 0, -P<Q<P\ (61) 

"Q > 0, -P < e < V , (62) 



The application of the formulas @ and Q together with the following equalities for the time derivatives 
of the Hubble parameter ( 15 1> 

H = -H 2 {q+l) , (63) 

H = H 3 {j + 3q + 2) , (64) 

H = -H A [k + 4j + 3q {q + 4) + 6] , (65) 

H = H 5 [l + 5k + 10j{q + 2) + 30q{q + 2) + 24} , (66) 
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allows us to write the higher-order energy conditions in terms of statefinders. The first-order dominant 
energy condition now reads as 



g + p = -2H A 



3q + 2 



K 



g-p = -2H* -j + 3g + 4 + 4 



a 2 H 2 
K 



~a 2 H 2 



the second-order reads as 



g + P 



2H 4 



k + 4j + 3g(<7 + 4) + 6 + 2- 



>0, 

>o, 

K 



(9 + 3) 



>0, 



g-p = 2H 

while the third-order reads as 

"Q + P = 



k + 2j + 3q(q + 6) + 12 + 4 



K 
a 2 H 2 



(9 + 3) 



> 



(67) 
(68) 

(69) 
(70) 

(71) 



g-p = 



2H b 



2H b 



1 + 5fc + 10 j(q + 2) + 30q(q + 2) + 24 + 2 



K 
a 2 H 2 



(j + 9q+ 12) 



> 



-1 + k + 2](4q + 11) + 6{7q 2 + 17 q + 8) + 4 (j + 9q+ 12) 



K 
a 2 H 2 



(72) 

>0, 



and so on. From these relations one can easily see that it is not possible to fulfill any generalized dominant 
energy condition if any of statefinders j, k, I, etc. is singular. This is due to the fact that the signs in the 
appropriate expressions in front of statefinders are the opposite. This gives a conclusion that the violation 
of the higher-order dominant energy conditions can be good signals for the emergence of the generalized 
sudden future singularities. 

In a similar manner one is able to consider the higher-order null, weak and strong energy conditions as 
follows 



g+p = -2H 3 
g+p = -2H 3 



J + 3<z + 2 
J + 3q + 2 



K 
~d 2 H 2 

K 
'~d 2 H 2 



>0, 

> , g> 



3 P = -6H d (j + q)>0 , g + p= -2H A [ j + 3g + 2 + 2 
■•• = >o, 



K 



2 H 2 



(73) 
(74) 

> , (75) 



and so on. 

After having a closer look at them it emerges that these higher-order energy conditions are not always 
very useful in order to determine the future fate of the Universe. 

However, some hybrid energy conditions may work for that. For example an energy condition 1 26 1 



ag > p, 

with a = const., which in terms of statefinders gives 



J > 



1 - 



3a \ 



2H 



: 1 



K 

' 2 H 2 



(76) 



(77) 



may prevent the emergence of sudden future singularity singularity for iV = 2 in 1441 . 
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6 Conclusions 

In view of the discussion performed in this paper one has the following remarks. 

Firstly, the future state of the universe may be more sudden and violent. It means that the universe may 
terminate either in a sudden future singularity or in a Big-Rip singularity and this is totally different from 
our earlier expectations that it could terminate in an asymptotically de Sitter state or in a Big-Crunch. 

Secondly, these new future singularities (Big-Rip and SFS) should not be confused - they have totally 
different properties with respect to geodesic completness. In particular, one can extend the evolution of the 
universe through a sudden future singularity, while it is not possible to do so for a Big-Rip. 

Thirdly, statefinders (Hubble, deceleration, jerk, kerk etc.) may be useful to diagnose the future state of 
the universe. By this we mean the emergence of sudden future singularities, the emergence of generalized 
sudden future singularities, the evolution of the cosmic equation of state w = w(t) etc. 

Finally, the new energy conditions may be introduced for the sake of the proper signal for generalized 
sudden future singularities, or (on the contrary), for the sake of the avoidance of sudden future singularities, 
or generalized sudden future singularities. 
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